Theoretical Study of Decomposition Pathways for HArF and HKrF 
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Abstract 

To provide theoretical insights into the stability and dynamics of the new rare gas 
compounds HArF and HKrF, reaction paths for decomposition processes HRgF -^Rg + 
HF and HRgF — »H + Rg + F (Rg = Ar, Kr) are calculated using ab initio electronic structure 
methods. The bending channels, HRgF -^Rg + HF, are described by single-configurational 
MP2 and CCSD(T) electronic structure methods, while the linear decomposition paths, 
HRgF — >H + Rg + F, require the use of multi-configurational wave functions that include 

dynamic correlation and are size extensive. HArF and HKrF molecules are found to be 
energetically stable with respect to atomic dissociation products (H + Rg + F) and 
separated by substantial energy barriers from Rg + HF products, which ensure their kinetic 
stability. The results are compatible with experimental data on these systems. 
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Introduction 

Very recently, a first stable chemically-bound compound of argon, HArF, was 
prepared experimentally in a low-temperature Ar matrix 1 ’". The corresponding krypton 
compound, HKrF, is also new and its preparation has just been reported 3 . Several 
theoretical studies have been performed that addressed the electronic structure and the 
nature of bonding in these molecules 4 ’ 7 , as well as their vibrational spectroscopy both in gas 
phase 5,7 and in rare-gas matrix environments 8 . Transition states for HRgF -^HF + Rg (Rg - 
Ar, Kr) reactions have been obtained previously 4 . In addition, stability of HArF with 
respect to the two dissociation channels (HArF — >H + Ar + F and HArF — >HF + Ar) has 
been considered by studying approximate minimum-energy reaction paths and energetics of 
dissociation products 6,7 . However, true minimum-energy paths for these processes have 
never been calculated. This study intends to fill this gap by performing intrinsic reaction 
coordinate (IRC) calculations following the minimum energy paths (MEPs) that connect 
HArF and HKrF compounds to their dissociation products. Knowledge of reaction paths is 
very important for understanding formation and decomposition dynamics of these novel 
rare-gas compounds. A similar study has been recently performed for HHeF compound 9 
that showed the importance of calculating reaction paths using the IRC method. By 
theoretical predictions, the species HHeF is a (first) chemically-bound compound of He in 
the electronic ground state 4 ' 5,9 though this is not yet experimentally confirmed. HHeF was 
found to be a metastable species separated from He + HF and H + He + F dissociation 


3 


4 


products by low (8-10 kcal/mol) barriers at the level of second order perturbation theory. 
Tunneling rate to decomposition of HHeF through these barriers has been estimated, and a 
relatively long lifetime (on the order of 120 ps) to He + HF decomposition has been 
predicted 9 . A study by Takayanagi and Wada at the multi-reference Cl level 10 showed lower 
(about 5 kcal/mol) barrier height for HHeF — »H + He + F reaction and shorter lifetime for 
HHeF compound (only 157 fs). This shows a great sensitivity of the predicted dynamics 
and behavior of the rare-gas compounds on the calculated potential energy surface. 

Although HArF and HKrF are more stable than HHeF, reliable estimation of the pathways 
between minima and dissociation products is crucial for understanding both dissociation and 
formation dynamics for these compounds. The calculated reaction paths will help future 
studies of the influence of temperature and IR excitation on decomposition dynamics and 
tunneling lifetimes of HArF and HKrF compounds. 


Computational methods 

Equilibrium structures of HArF and HKrF minima and of transition states (TS) 
leading to HF + Ar(Kr) products were computed at the second-order Moller-Plesset (MP2) 
level of electronic structure theory 11 with augmented correlation-consistent valence triple-^ 
(aug-cc-pVTZ) basis set 12 ' 14 . The same level of theory was used to obtain the minimum 
energy paths (MEPs) leading from the equilibrium structures to the energetically more 
stable products Ar(Kr) + HF. These MEPs were calculated using the Intrinsic Reaction 
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Coordinate method (IRC) 15 with the second order Gonzalez-Schlegel algorithm and a step 

size of 0.1 amu 1/2 *bohr. Single point energies were estimated along these MEPs at the 
higher level coupled-cluster method CCSD(T) 17 . Stationary points optimizations and IRC 
calculations were carried out using the electronic structure package GAMES S , while 
CCSD(T) single-point energies were computed using GAUSSIAN98 19 . 

The IRC paths to H + Rg + F were studied using CASSCF wave functions that 
included 16 active electrons and 9 active orbitals. This is the full-valence active space 
(excluding valence d-orbitals) and is denoted as MCSCF(16,9). The same, aug-cc-pVTZ, 
basis set was used as for the bending (HF + Rg) channel. The potential energy points along 
the IRC were computed using several methods that include dynamic electron correlation 
effects: multi-reference second order perturbation theory CASPT2 20 , internally contracted 
multi-reference configuration interaction MRCP 1 , and size-extensive modifications of 
MRCI: multi-reference averaged coupled-pair functional MR-ACPF 22 and multi-reference 
averaged quadratic coupled-cluster AQCC" 3 methods. CASPT2, MRCI, MR-ACPF, and 
MR-AQCC calculations are all based on MCSCF(16,9) reference wave functions and 
performed using MOLPRO electronic structure package 24 . 

Results and discussion 

The reaction paths for the HArF -^>Ar + HF and HKrF -*Kr + HF dissociation 
processes are shown in Figs. 1 and 2. The IRC paths were obtained at the single - 
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configurational MP2/aug-cc-pVTZ level of theory. The validity of the single-configurational 
approach for this dissociation channel was tested in our previous study for a similar 
reaction path for HHeF 9 . It was shown that multi-configurational second order perturbation 
theory (CASPT2) gave virtually the same potential values along the reaction path as smgle- 
configurational MP2. Therefore, we assume that single-configurational MP2 approach 
should be valid in the case of HRgF — >Rg 4 HF dissociation channels also for Ar and Kr. To 

test the MP2 reaction path further, we have calculated CCSD(T) potential energy points 
along the MP2 reaction paths. It can be seen from Figures 1 and 2 that CCSD(T) relative 
energies along the IRC paths are in very good agreement with and practically 
indistinguishable from MP2 results. This also supports the validity of the MP2 approach 
for HRgF — >Rg + HF dissociation channels. The transition states found in this study are 

23.7 and 32.1 kcal/mol above the equilibrium structures of HArF and HKrF, respectively. 
These are slightly lower than predicted by Wong (Ref. 4): 28.0 and 34.6 kcal/mol. The 
configurations corresponding to the transition states are shown in Fig. 3. The IRC paths in 
the beginning correspond approximately to the motion along the bending coordinate of 
HArF (HKrF). At later stages, the rare gas atom assumes a closed shell electronic structure, 
and the H-F bond forms (see Fig. 1). The dissociation products (Rg + HF) are lower in 
energy than HRgF by 134 kcal/mol and 1 12 kcal/mol in case of Ar and Kr respectively. 
However, relatively high barriers to these dissociation products make both HArF and HKrF 
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compounds kinetically stable. This is compatible with the experimental preparation of both 
HArF and HKrF 1 ' 3 . 

The reaction paths for the 3-body decomposition channels HArF -> H + Ar + F and 
HKrF — > H + Kr + F are shown in Figs. 4 and 5. The dissociation motion in this case is 

collinear rather than bending. It can be seen from these figures that MCSCF method 
incorrectly predicts both HArF and HKrF to be energetically higher than atomic 
dissociation products. Obviously, more accurate methods that include dynamic correlation 
are necessary to obtain more reliable results for these dissociation channels. Although the 
situation is greatly improved at the MRCI level as compared to MCSCF, this level still 
does not provide quantitatively correct 3-body dissociation paths and predicts dissociation 
limits too low in energy with respect to HRgF minima. CASPT2 calculations predict both 
dissociation products and transition states at higher energies with respect to HRgF minima 
than MRCI. At the level of CASPT2, the dissociation products are predicted to be about 
6.8 and 26.4 kcal/mol above the HArF and HKrF minima respectively. This quite significant 
difference between MRCI and CASPT2 results has prompted us to test additional methods 
that combine advantages of both MRCI and CASPT2: size extensive multi-reference ACPF 
and AQCC methods. As can be seen from Figures 4 and 5, these methods give energies, 
which are in between those of MRCI and CASPT2. They predict both HArF and HKrF 
compounds below their dissociation products H + Rg + F, which is an encouraging result. 
The binding energies for HArF are 0.2 and 1.8 kcal/mol from AQCC and ACPF 
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respectively. This is consistent with experimental results 1 ' 2 that show that HArF is stable 
with respect to H + Ar + F, but the binding energy is not very large. The binding energies 
for HKrF are 19.8 and 21.2 kcal/mol at the AQCC and ACPF levels respectively. These are 
somewhat lower, but not very different from the CASPT2 result. It should be noted that 
there is a sizable barrier for the formation of HArF from atomic H + Ar + F fragments, 
about 10-11 kcal/mol. This barrier appears to be too high since the HArF can be formed 
upon annealing of solid argon matrix at about 1 8 K 1 . On the other hand, the large bamer 
from HArF towards its collinear dissociation channel seems reasonable in the light of 
annealing experiments of HArF in solid Ar 2 , in which the IR spectrum of HArF in a 
thermally relaxed environment is visible until the sample evaporates above 40 K. In the case 
of HKrF the dissociation barrier of the collinear H + Kr + F -> HKrF channel is much 
lower, about 5-6 kcal/mol. This correlates well with the experimental observation 3 of HKrF 
formation in one of the sites of Kr matrix at about 13 K, which suggests a reaction with a 
very low barrier. It can be concluded that potential energy surfaces for HRgF — >H + Rg + F 
dissociation channel are not easy to describe correctly for such unusual rare-gas compounds 
as HArF and HKrF. MCSCF and internally contracted MRCI methods do not give 
qualitatively correct energetics for this reaction channel. On the other hand, CASPT2, MR- 
AQCC, and MR-ACPF calculations seem to reproduce the experimentally observed 
features of formation of these two molecules qualitatively correctly and could be considered 
a rather good starting point for understanding both dissociation and formation dynamics on 
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molecular level in solid hosts. Work on tunneling decomposition and on formation dynamics 
of HArF, HKrF in the corresponding Ar and Kr matrices is currently under way 

Conclusions 

Minimum-energy reaction paths are computed for decomposition processes HRgF 
^Rg + HF and HRgF — >H + Rg + F (Rg = Ar, Kr) using ab initio electronic structure 

potentials and intrinsic reaction coordinate (IRC) method. Although some of these reactions 
have been studied previously, true minimum energy paths have not been obtained before. 
The calculated energetics for HRgF -^Rg + HF reactions (barrier heights and 
exothermicities) are in good agreement with previous studies. While the HRgF ->Rg + HF 
channels are described reasonably well by single-configurational MP2 and CCSD(T) 
electronic structure methods, investigation of the linear decomposition paths (HRgF ->H + 
Rg + F) is more challenging and requires use of multi-configurational wave functions, as well 
as dynamic correlation treatment. Furthermore, different methods that account for dynamic 
correlation (CASPT2 and MRCI) give quite different results for this reaction path. It 
appears that the results obtained using size extensive multi-reference ACPF and AQCC 
methods provide the most reliable predictions in this case. The shapes of the reaction paths 
for HArF and HKrF dissociation processes are obtained in this study for the first time and 
should be useful for future studies of the influence of temperature and IR excitation on 
decomposition dynamics and tunneling lifetimes of HArF and HKrF compounds. 
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Figure Captions 

Figure 1 . Minimum energy path for HArF — > Ar + HF reaction. 

Figure 2. Minimum energy path for HKrF — > Kr + HF reaction. 

Figure 3. Equilibrium and transition state geometries for HArF and HKrF obtained at the 
MP2/aug-cc-pVTZ level of theory (bond distances are given in angstroms, angles - in 

degrees). 

Figure 4. Minimum energy path for HArF — » H + Ar + F reaction. 

Figure 5. Minimum energy path for HKrF — » H + Kr + F reaction. 
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